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Tj- ' Abstract 

The fluid shear stress is an important regulator of the cardiovascular system via the endothelium- 
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fH , 

O . derived relaxing factor (EDRF) that is Nitric Oxide. This mechanism involves biochemical reactions 

in an arterial wall. The autoregulation process is managed by the vascular tonus and gives the 
negative feedback for the shear stress changing. A new mathematical model for the autoregulation 
of a blood flow through arteria under the constant transmural pressure is presented. Endothelium- 

(— I , derived relaxing factor Nitric Oxide, the multi-layer structure of an arterial wall, and kinetic- 
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diffusion processes are taken into consideration. The limit case of the thin-wall artery is analytically 
studied. The stability condition for a stationary point of the linearized system is given. The exact 



S2^ . stationary solutions of the origin system are found. The numerical simulation for the autoregulation 

.^ I system is presented. It is shown the arteria adaptation to an initial radial perturbation and the 

O ' 

\^ ' transition of the system to new equilibrium state in response on the blood flow changing. 
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I. INTRODUCTION 

The modeling of a cardiovascular system is the problem of great importance due to 
it should be aid for understanding and prediction of various diseases like atherosclerosis, 
arteriosclerosis, hypertension, etc. The most often approaches to the problem is a direct 
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^ plying of classical hydrodynamics models of fluid flow through elastic shells or tubes 
3|. But in some cases, for example, for muscular resistance arteries, it is necessary to take 
into account the difference between an ordinary passive and a "biological" active tube. 

The one way to reach the purpose is a wide discussed effect of the flow-induced vasodila- 
tion by Nitric Oxide (NO) radical |^, |5|, la l7[ . For a long time the endothelium cells, covering 
the arterial bed surface, was supposed to provide only the friction reduction for the blood 
flow through arteria. But since the EDRF was discovered via the comparison of two arterial 
rings with and without endothelium by the ability to acetylcholine-dependent smooth mus- 
cle relaxation and due to later investigations it was shown the endothelium plays the main 
role for the local blood flow regulation. The so-called Endothelium Derived Relaxing Factor 
(EDRF) was discovered in 1980 by Robert F. Furchgott [8|, |9|. Nitric oxide was proposed 
as a signal molecule to set a connection from endothelium to the smooth muscles. The 
EDRF-NO mechanism enabled to explain the principle of action for the flrst-aid medicine 
Nitroglycerine being used before without understanding. 

The mechanical nature of an arterial wall tonus regulation is highlighted in the works 



12| |. It is perceived that the increasing of a shear stress between the blood flow and 
inner arterial surface causes the relaxation in a smooth muscle layer of an arterial wall. This 
necessary induces the increasing of the arterial radius and the decreasing of the shear stress 
itself. Therefore the process in a whole gives the system with a negative feedback. 

There are three main layers in an arterial wall. The flrst, internal layer is intima (i), 
the second layer is media (m) and the last one is adventitia (a). The inner boundary of 
the intima layer, i.e. the internal arterial surface, is covered with endothelium cells. The 
media layer is full of smooth muscle cells. The thicknesses of the layers depend on the type 
of artery or arteriole. We mainly consider the muscle resistance arteries which have well 
developed muscle layer (media) and non-vanishing intima layer. The typical ratio of intima 
thickness to the media one is about 10^^. 

The scenario of the flow-induced relaxation is as follows. The increasing shear stress cr shear 



on the surface of endothelium cells opens calcium channels which launch the production of 
NO from L-arginine with NO-synthase (NOS) catalysis and then NO diffuses with descending 
through the intima layer towards the smooth muscle cells in the media layer. As a lipophilic 
molecule NO easily penetrates through a cell membrane of a muscle cell and initiate synthesis 
of the cyclic guanosine monophosphate (cGMP). Ultimately, cGMP stimulates the outflow 
of intracellular Ca^^ that leads to relaxation of the smooth muscle cell. The flow-induced 
contraction is realized vise versa. 

The aim of this paper is to develop and to study the mathematical model for description 
of blood flow autoregulation that accumulates a viscoelastic nature of an arterial wall and 
the two-layer diffusion and kinetic processes for concentrations of the key agents: Nitric 
Oxide (NO) and Calcium ions {Ca^^). 

The outline of the article is as follows. In the section |n] we introduce the assumptions of 
the model. In the section ITTTl we derive the closed system for description of the autoregulation 
process. In the section IIVI the steady-state concentrations of NO and Ca?^ are obtained. 
In the section El we study the limit case of a thin-wall artery. In this case the stability 
condition of an equilibrium state of the system is given. In the section IVII we consider the 
case of passive dilation of an arteria with the fully relaxed muscles. The exact kink-shaped 
solution is found. In the section I VI II the numerical simulation of the autoregulation process 
near the stationary state is presented. In the section IVIIII we summarize and discuss the 
obtained results. The appendix El includes the essential notations. The appendix JBl gives 
the approach for finding exact solution of the passive vessel model. 

II. MAIN ASSUMPTIONS OF THE MODEL 

We consider the arteria to be axial-symmetric, viscoelastic and incompressible. The 
blood is also assumed to be incompressible and Newtonian. The flow is quasi-stationary, 
the transmural pressure is constant and the velocity profile is the power generalization of 
the Poiseuille's law. We suppose the dependence of muscular force on calcium concentration 
to be linear and the dependence of Ca^^ concentration decreasing ratio on Nitric Oxide 
concentration is also linear. The concentration of NO in the endothelium is assumed to be 
proportional to the shear stress on an arterial wall jSi]- 



III. THE STATEMENT OF THE PROBLEM 

Let us consider an arterial segment of length / in the cylindrical coordinate system 
(r,9,x = z). The intima, media, and adventitia layers have coordinates i?j, /?„, Ra re- 
spectively. 

A. The shear stress dependence on the blood flow 

Consider the power generalization velocity profile of the Poiseuille's flow J3|: 
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Here Vx is the axial velocity, u is the cross-sectional averaged axial velocity, R is the arterial 
radius and 7 is the profile sharpness. 

In case of Newtonian fluid with dynamical viscosity /i the shear stress on the wall of 
elastic tube is 
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where Q = Au is the blood discharge through the cross-section with the area A. 

Summarize all the assumptions for the laminar stationary flow, the cross-section averaged 
Navier-Stokes equation takes the form of the generalized Hagen-Poiseuille equation: 

^^2h.2„-|, (3) 

where AP is the pressure difference on an arterial segment with the length /. 

It shows the linear dependence of the pressure gradient from the discharge and inversely 
proportionality to the fourth power of the arterial radius. 

In case of axial-symmetric radial perturbations R{t) = Rq + r](t) we have from (j21): 
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There is a hypothesis of maintaining the shear stress principle ashear = const J6|,|l3|. One 
can conclude the increasing of the flow necessitate the increasing of the steady-state arterial 
radius to compensate the changing of the shear stress. The estimated relation between the 
new steady-state discharge and new stationary radius is as follows: 




^=(7^-1^0 (5) 



It is remarkable the difference of reaction for increasing and decreasing the blood flow 
near the previous stationary value. The changing of the radius in response to higher flow is 
smaller than for the same lower flow. It is explained by the inversely cubic dependence of 
shear stress from the radius. 

One can see there is the linear dependence between radial perturbation and mean blood 
flow in case of small radial perturbations (|?7| << Rq). 

B. The synthesis and diffusion of Nitric Oxide 

According to the EDRF-mechanism mediated by the fluid flow, the concentration of Nitric 
Oxide produced by an endothelium cell is managed by the shear stress value. We consider 
the NO transport to the smooth muscle tissue as a diffusion process (diffusion coefficient is 
Di) with a descending (reaction coefficient is 6i). Then it continues to diffuse through the 
media layer but with another diffusion coefficient D2 and reaction coefficient 62- 

The production of Nitric Oxide in an endothelium cell has the shear stress (Jshear as one 
of essential regulators, therefore this process can be described with a kinetic equation: 

-keUe + lpashearit) (6) 



dt 

where Ue is the NO concentration in an endothelium cell, kg is the rate of mass transfer of 
NO from the cell, and ip is the production rate constant. 

Under the assumption of quasi-stationary NO production, i.e. that characteristic time 
of NO ma=s truster from an endothelium cell towards intima law (r„o_.„„„_.„.„,„ ~ 

Ar'^/D ~ 1/528 sec, where D = 3300 jjm'^ / sec, Ar ~ 2.5 /im |18|) is smaller than the 
typical time of ashear changing {Tshear ~ Tradius-osdUations ~ 1/2 scc), from equation © we 
have the following relation between the rig and cr shear'- 

ip 

IT-eit) = TT^^shearit) (7) 

The relation ((7j) is used as inner boundary condition for Nitric Oxide diffusion through an 
arterial wall (^|T-=_Ri„ti,„a = ^e)- Ultimately, at the inner boundary of intima layer we assume 
the concentration of NO to be proportional to the shear stress (proportionality coefficient 
is fcs). Between the intima and media layers we use the continuity of concentrations and 
fluxes. On the external layer we take the impenetrability condition into account. Thus the 
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system of equations and the boundary conditions for the Nitric Oxide concentration are as 

follows: 

^ = D-^(r^]-6n-, 
dt ■' r dr \ dr J ■^ "" 

Ri < r < Rm for j = 1 (intima) 
Rm < r < Ra for j = 2 (media) 

ni\r=Ri = kz Cr shear 



r=R„ 



The system of equations (jH)) together with initial conditions describes the two-layer 
diffusion-kinetic process for the Nitric Oxide in an arterial wall. 

C. The equation for the kinetics of the Calcium ions in a smooth muscle cell 

To derive the balance-equation for concentration of Ca^^ in a smooth muscle cell it 
is necessary to describe the ways of Ca^'^ in- and out-fiuxes. There are two source of 
the calcium ions: the extracellular space and the intracellular containers - sarcoplasmic 
reticulum. The concentration of Ca^^ in these sources is about 10^ greater than in the 
intracellular space. The balance of calcium ions in the muscle cell at the point r may be 
described, similarly to [gI], as 

^^Tl^ = -a{C - Co) + Pia^t -C)~ hn2ir, t) (9) 

ot 

where the first term is responsible for the natural active outflow transport of Ca"^^ com- 
pared to the minimal observed concentration Co, the second term is described a passive 
diffusion provided by the difference between the intracellular calcium concentration C and 
extracellular ones Cext, and the last term is presented the NO-mediated active outflow. 

Taking into account the relation Ce^t » C we can treat it as a constant source: ipo = 
aCo + pCext = const. In this case the equation Q can be transformed to the form: 

^^^ = -aC - hn2 (r, t) + ^o (10) 

The equation pO|) is used to describe the Calcium-balance in the smooth muscle layer. 



D. The equation for an arterial wall movement 

In order to obtain the close system of a blood flow autoregulation we need to have a 
link between the radial perturbation and the external forces such as pressure and muscular 
force jl^. The constitutive equation 15| can be found from the movement equation for an 
arterial wall segment. 

Let us consider the incompressible viscoelastic wall element with mass Am, density p^, 
width h, radius R, and length Ax. According to the movement law 

A ^- f ^f 
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dt^ (11) 

f radial = -(^09 27TAxh, f pressure = (-P " Pext) ^TlAxh 

where Am = p„,27ri?Ax/i, f radial is proportional to the circumference component of a stress 
tensor age and /pressure is the resulting transmural pressure (the difference between the 
internal and external pressure). 

The stress tensor component aeg consists of three parts: a passive elastic force (weakly 
nonlinear with quadratic addition), a viscous force and an active force due to the muscle 
tonus 

--^l^^-^'(^)^-^f-^-=^(^' (X2) 

here E{< C >) is the Yung's modulus dependent on averaged concentration of Co?'^ in a 
muscle cell layer, ^ is the Poisson's ratio, Ei is the small nonlinear elastic coefficient for a 
square addition, A is the viscous characteristic of the wall, F{C) is the active force component 
determined also by the integral calcium concentration level above the threshold one Cth-, ^2 
is the coefficient of proportionality for the muscular tonus response on the Ca?^ level. 

Substitute (fT^ in (fTT|) and take into consideration the linear dependence of muscle force 
on calcium and the incompressibility condition /lo-Ro = hR. Then the constitutive equation 
for the radial perturbations {R = Rq + r], \ri\ << Rq) has a form 
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where 

x(C)=xo(H-eF(C)), xo= f° £«1 

^o(l - 4 ) 

F(C^)= /" \c-CtH]e{C-Cth)rdr, (14) 

^- the Heaviside's step function 

Renormahze the constants A, x, k2 with the value ho/Ro and denote the constant, under 
the assumptions, transmural pressure Pq = P — P^xt = const and Xi = -^l^. 

Ultimately, we obtain a new integro-differential equation describing the wall movement 
in the presence of smooth muscle tonus 

P«,/io^ + A^ + x(C)77 + xi772 = Po- k2 I \c-CtH]e{C-Cth)rdr (15) 

One can see in the case of absence of muscle force (full relaxation) it is the equation of 
a nonlinear damping oscillator with an external force. The presence of calcium-dependent 
force term is provided the feedback and makes the arteria different from a passive viscoelastic 
tube. 



E. The problem statement for the blood flow autoregulation in dimensionless 
variables 

Summarize the equations obtained above we have the complete system to describe the 
process of blood flow autoregulation due to EDRF-NO mechanism: 
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with the boundary conditions: 
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As the initial values it is taken the perturbed steady-state solutions. 

Here equation (fTT)|l describes the Ca^"'"-balance in a smooth-muscle cell, equations 
(fTTj) . (fTHj) characterize the diffusion of Nitric Oxide in intima and media respectively, and the 
equation ()19|) gives the relation establishing the arterial wall movement under the influence 
of the average calcium ions concentration. 

In order to pass on to the non-dimensional system of equation setting up the new dimen- 

sionless variables: 

C = Cth C, ni = n? n[, n2 = n" n'^, 

r] = r]or]', t = to t', r = r^r' 
were for convenience choosing 
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After substitution (PT|) the system (fTBj) - ()19|) turns into a dimensionless form (primes 
over the variables are omitted): 
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where dimensionless constants are 
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Then the boundary conditions take the form: 
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where i?o = Rti R'm = Rm/Ro, R'a = Ra/Ro and initial values are the perturbed 
solutions of the steady-state system. 

From the non-dimensional system of equations one can make a remark that the stationary 
blood flow discharge through the vessel's cross-section Q has implicit influence on the Ca'^'^- 
concentration in the smooth muscle cell via term k[ n2 in the equation (j^H|l due to the 
coefficient k[ r^ n^ ^ Q. 

IV. THE SOLUTION OF THE PROBLEM IN A STEADY STATE 



To consider the stationary case letting the following: 

C = C{r), rii = hi{r), n2 = n2{r), R = Ro = const 
Under the assumptions the system of equations (j23p - (j26|) takes the form: 

C(r) = -k[n2ir)+^o, R'm < r < K 



(29) 
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with the boundary conditions: 
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The ODEs (j3ip . (j32j) for NO concentration have general solution via modified Bessel 
functions lQ{z),KQiz): 
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where Ai, A2,Bi, B2 are the arbitrary constants defining by the boundary conditions ()34|) : 

^1/0(6) +A2iro(6) = l 

5i/i(6K)-52i^i(6K) = 
A,Ui,IiJ + A2Ko{iiE!J = 

i,{A^h{^,B!J-A2K,{i^IiJ) = 

where ^i = ^J§r^, 6 = ^|%f^ 

Using the typical experimental data for the muscular resistance artery |l6l . llTl |l8| 
Ri = 1.0 mm, h = 0.5 mm; R'^ = 1.05, R'^ = 1.3 and assuming ^1 = 6, ^2 = 2 we can find 
the constants Ai, A2, -Bi, -B2 from the boundary conditions ()36|) . 

Steady-state Ca'^'^- concentration C{r) is given by (jHUjl . 

The equilibrium distribution of concentrations is depicted on the figure d 
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FIG. 1: Stationary distribution of NO and intracellular Calcium ions. 



V. THE CASE OF A THIN- WALL ARTERY 



To understand the qualitative behavior of the system consider the limit case of a thin wall 
artery. The similar model was studied by A. Rachev, S.A. Regirer et al. in J6|,ll9|. There 
are estimations to come to the limit case. The first relation is hi/hm « 1 that enables to 
come to a one-layer wall model. The second one is r^if fusion « Tkinetic, were the typical 



time of the diffusion process is r^. 



'if fusion 



jT and the typical time of the kinetic process is 



D 



Tkinetic = 'min{^ , ^}. Here hi, h^ are the wall thickness of intinia and media layers, h is the 
spatial scale of the wall thickness. Considering that the kinetic processes for Nitric Oxide 
are faster than for Calcium ions we have as follows: 



h « 



^ -h 



(37) 



where ho is the characteristic wall thickness to compare with. Taking into account the typical 
values of the parameters as Z^ = 3300 yum^/sec and 6 = 1 sec~^ [18] we obtain h^ = 57 fim. 

It is also should be note the default condition of quasi-stationary diffusion: T^if fusion « 
Tosc, where Tosc is the typical period of radial oscillations. The typical value of T~l is about 
1 -h 2 sec~^ then the /iq value is close to 57 fim or a bit less. 

The large and medium resistance muscle arteries have the specific wall thickness 
h ~ 100 -7- 1000 fim whereas the small arteries and arterioles have the much smaller thick- 
ness /i ~ 10 ^m. Therefore the limit case covers the case of flow in a small artery with 
h « 50 /xm. 
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Thus the intiina and media layers are so thin to neglect the multi-layer nature of the wall 
and eliminate the diffusion processes. 

After the averaging of the calcium feedback F{C) over the wall thickness the system 
(I16|) - (jl9|) takes a simplified form: 

dx a 

— — = —a X n + 

dt (1 + 1)3 

dy ^ ^ (38) 

dt 

dz r> 2 7-> 

— — = —A X — ny — pz — Kiy — K2xy + B 
dt 

were x = x{t) = C{x,t) — Cth is the average concentration of Ca^^ in the arterial smooth 
muscle layer, y = y{t) = rjit) is the deviation of the radius of the vessel {y > — c), z = z{t) 
is the velocity of radius oscillation; a is the rate of a natural "pumping" of the free calcium 
ions from the intracellular space, a is represents the blood flow level (a ~ Q), c is the 
non-perturbed arterial radius, b is the rate of the calcium inflow in a smooth muscle cell, A 
is the coefficient of proportionality for the calcium-feedback force, k is the linear elasticity 
coefficient, ki is the nonlinear elasticity coefficient, k,2 is the small calcium- induced elasticity 
coefficient, (3 is the viscous (resistance) coefficient of an arterial wall, B represents the mean 
constant transmural pressure. 

Look for the stationary points of the system ()38|) . One can obtain under the condition 



B , , 

b — a = a— (39) 

there is a stationary point {x = B/A, y = 0, z = 0}. It corresponds to the non-perturbed 
state of an artery. All the rest real stationary points of the system have y < —c and 
hence they are out of physical sense. The relation (jHIJjl reflects the balance between the 
muscle forces mediated by calcium concentration and the pressure forces in the blood. The 
steady-state Ca^'^ concentration is equal to x = B/A ~ Po/(^o-Ro)- 

Study the stability of the dynamical system (jHHI) near the stationary point {B/A, 0, 0} 
taking into account relation ()39p. Consider the linearized system 

§ = AA' + F, (40) 
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X = {x,y,zf, F={b-a,0,Bf 
The Routh-Hurwitz criterion provides the condition then all eigenvalues of A have nega- 
tive real parts. Here the stability condition is as follows: 



P ( a^ + a/3 + K + K2— ) > 



3aA 



(41) 



Taking into consideration the strictly positiveness of the A, a, c, a and non-negative values 
of the rest parameters one can conclude from (PT|l the condition for the wall viscosity (3 > 0. 
It shows the importance of the viscoelastic nature of an arterial wall to maintain the stability 
of the stationary state. In general case, there is the critical wall viscosity PcriUcai below that 
oscillations demonstrate the lack of stability. 

The qualitative analysis on the phase plane confirms the preliminary estimates (figure |2I). 




FIG. 2: The two-dimensional projection of the phase trajectory of the system. For (3 = PcriUcai it 
is the periodic oscillations (left) and for /3 > ^critical it is the damping oscillations (right). 



VI. THE CASE OF A PASSIVE VESSEL 

One can see from the thin-wall approximation the more the discharge the less the equi- 
librium calcium level. In the general model we have a non-constant distribution of Ca^^ 
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concentration. If the stationary Ca^^ concentration for the whole arterial wall is below the 
threshold level Cth it becomes a fully relaxed. In this case the 'active' viscoelastic tube is 
reduced to the 'passive' one. The law of the arterial wall motion (j26j) in the dimensionless 
form (primes are omitted) is as follows: 

d Ti d Ti 

_| + A-^ + x„„ + x,„^ = P„ (42) 

The nonlinear differential equation (|42p can be solved exactly via the simplest equation 
method J2iJ, |2l| . One can obtain 

r,{t) = f]^ tanh (^] (2 - tanh (^) ) (43) 



Po /4xiPo , 4/2500 xiPo 



3xi V 3 V 27 

The kink-shaped solution demonstrates the switch from one steady state to another under 
a constant force field. 

Here the solution (jIHj) satisfies a non-perturbed state of artery with rj{Q) = 0. The 
pressure and smooth-muscle force compensate each other. After a vanishing of the muscle 
force (due to the sharp decreasing of the calcium level) arteria expands to a new equilibrium 
state. 

The new arterial radius depends on the transmural pressure and the elastic properties of 
an arterial wall. It can be estimated by r]oo- 

VII. THE NUMERICAL SIMULATION FOR THE PROBLEM OF BLOOD FLOW 
AUTOREGULATION 

Consider the general case of the two-layer kinetic-diffusion system in the dimensionless 
form (j23p - (|28|) for description of the blood flow regulation. In order to study the dynamics 
of the solutions of the system near the steady state the numerical simulation is performed. An 
implicit iterative finite-difference scheme is implemented. As the initial values the perturbed 
exact stationary solutions (IHUj). (jHSj) are taken. 

The behavior of the solution for the initial stretching of the radius confirms the asymptotic 
stability of the stationary state (figure E))- 

As a test solution in case of passive dilation the exact solution ()43j) is taken. 
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FIG. 3: The dynamics of the system relaxation to the previous steady state after an initial stretching 
of the artery r]{t = 0) = 0.1. 

The comparison gives a good agreement between the numerical solution and the exact 
one (figure EI). 
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FIG. 4: The passive expanding of the arteria due to a constant transmural pressure. The exact 
solution (|15|l (dotted line) and the numerical one (solid line). 

In response on the changing of the blood flow, that is managed by the coefficient ki ~ Q, 
the system comes after the damping oscillations to a new steady state (figure Ej). It is 
remarkable the different reaction of the system to the increasing and decreasing of the 
discharge. The relaxation time in case of fiow decreasing is smaller than in case of fiow 
increasing. It may be explained by the drop of the critical viscosity level in response on 
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the decreasing flow according to ()4H) . Also the arterial radius deviation is bigger in case of 
decreasing of blood flow accordingly to the inverse cubic dependence of the shear stress on 
the radius. In case of increasing flow it is vice versa. 

One can see the growth of the blood flow can potentially be a source of instability espe- 
cially for small arterial wall viscosity near the critical one. 




90 100 



FIG. 5: The transition of the system to the new equilibrium state after the decreasing for 25% 
(left) and increasing for 25% (right) of the mean blood flow. 



VIII. CONCLUSION 



The two-layer diffusion-kinetic model is proposed to describe the process of a local blood 
flow regulation in an artery. The exact stationary distributions of the key agents - Nitric 
Oxide and Calcium ions are obtained. 

The limit case of a thin wall artery under the estimation of the wall thickness (I37|) 
is analytically studied. The stability condition for the equilibrium state is given by the 
formula (^1]) . It is shown the necessity of the viscoelastic nature (non-zero viscosity) of the 
arterial wall to provide the stability of the system. The minimal critical viscosity value of a 
wall is obtained in the linearized case. 

In case of full relaxation of the smooth muscles the exact solution in the kink form is 
found to describe the passive dilation of the artery. 

The numerical simulation demonstrates the transition of the system to the new steady 
state with the new radius value in response of changing of the mean blood discharge. This 
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result is in agreement with the experimental observation 4]. It is confirmed the importance 
of the endothelium derived relaxing factor - Nitric Oxide for arterial haemodynamics. 

The model can be applied to the study of the local autoregulation of the coronary, cerebral 
and kidney blood flow. 
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APPENDIX A: THE SUMMARY OF NOTATION 

Assuming the cylindrical coordinate system {r,6,x = z). 
Ri is the coordinate of intima layer boundary (equal to inner vessel radius); 
Rm is the coordinate of media layer boundary; 
Ra is the coordinate of adventitia layer boundary; 
R = R{t) is the radius of the arterial wall {R = Ri is the inner radius); 
Tj = rj{t) = R(t) — Rq is the perturbation of the steady-state arterial radius Rq; 
C = C{r, t) is the concentration of Ca^"*'-ions in the smooth muscle cell; 
Cth is threshold concentration of Ca'^^ to start the contraction in the smooth muscle cell; 
ni = Ui (r, t) is the concentration of A^O-radical in the first (intima) layer; 
n2 = n2{r,t) is the concentration of NO in the second (media) layer; 
Pq = P — Pext = const is the cross-section averaged stationary transmural pressure; 
u is the cross-sectional averaged axial fluid velocity of a steady-state flow; 
Q = Au = const is the fluid discharge through a cross-section of artery. 

APPENDIX B: FINDING EXACT SOLUTION OF THE NONLINEAR ODE (|42|l 

To obtain an exact solution of the equation: 

+ \-— + Koy + Kiy = Pq (Bl) 



dz"^ dz 



we use the simplest equation method [20| that genera 



izes the existing approaches like the 



tanh- method, the method of trial elliptic functions J2l|, etc. 
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Taking into account the second order pole of the general solution of ()B1|) look for solution 
in the form of the following expansion: 

y{z) = A^ + ArG{z) + A2G{zf (B2) 

where G{z) is the solution with the first order pole of the equation 

dG{z) 



dt 



kG{z)-kG{zy (B3) 



and Aq, Ai, A2, k are the arbitrary constants to be determined. 
Substituting the expansion ()B2|) into equation ()B1|) we find 



30AA; + 25x - A2 + 25A;2 _ Qk{X + 5k) 

Aq — , Ai — 

50 xi 5xi 

6A;2 ^ -36A^ + 625x2 , , A 

An = , Fr) = , k = ± — 

xi ' " 2500xi ' 5 



(B4) 



Taking the solution G{z) = | + |tanh(|/c(z — zq)) of the auxiliary equation ()B3|) and 
choosing /c = ^ we have the exact solution of the ODE (jBlj) in the form: 



y{z) = -^ (SX^ - 25x + GA^ tanh \^X(z - zq)] - 

50x1 V (25^ 

-3XHanh^[^X{z-zo)]\ 

where Zq is an arbitrary constant. 

Under the additional condition y{0) = 0, zq = there are the following relations between 
the parameters: 

After the simplification of ()B5|) taking into account ()B6|) we obtain finally the kink-shape 
solution of the equation ()B1|1 : 



viz) = y^ tanh ( ^) (2 - tanh (^^ ) ) (B7) 



where yoo = Jft^- 
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